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SECTION  I 

INTRODUCTION  AND  SUMMARY 


This  third  quarterly  report  summarizes  the  progress  made 
during  the  period  10  January  1976  to  10  April  1976  in  the  VELA  Network 
Evaluation  and  Automatic  Processing  Research  program  being  carried  out 
by  Texas  Instruments  Incorporated  at  the  Seismic  Data  Analysis  Center 
(SDAC)  in  Alexandria,  Virginia.  The  four  program  tasks  are: 

• Array  and  network  evaluation 

• Signal  detection  methods 

• Signal  estimation  techniques 

• Discrimination. 

During  the  last  quarter,  processing  and  analysis  of  long -period 
and  short-period  data  from  the  Albuquerque,  Guam,  and  Mashhad  SRO  stations 
was  initiated,  with  the  goal  of  establishing  the  noise  characteristics  and  signal 
detection  and  discrimination  capabilities  of  each.  The  quality  of  the  SRO  data 
sampled  has  been  in  general  excellent,  with  few  long-period  transients  and  no 
spikes  observed.  Software  modifications  were  completed  for  the  ILPA  eval- 
uation, and  this  evaluation  will  begin  when  the  data  are  available.  Work  has 
started  on  the  program  to  transfer  ILPA  data  to  the  mass  storage  center. 

Evaluation  of  the  automatic  power  and  Fisher  detectors  using 
KSRS  short-period  data  was  continued.  Detection  probabilities  using  37 
Eurasian  events  recorded  in  November  1974  were  obtained  for  several  pre- 
filters with  0.8  and  1.6  second  integration  gates.  For  optimum  detection  with 
the  Fisher  and  the  conventional  power  detectors,  a single  sensor  prefilter 
(1.  5 to  2.5  Hz)  has  been  determined.  During  the  next  quarter,  detection 
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probabilities  for  three  regions  will  be  examined  using  February  1976  (winter 
season)  data. 

A modified  version  of  the  program  NETWORTH  has  been  used 
to  evaluate  the  effect  of  certain  physical  parameters  on  network  detection 
capability.  This  study,  which  was  completed  during  the  last  quarter,  indi- 
cates that  the  station  detection  threshold  and  the  station  reliability  are  the 
primary  determinants  to  network  performance.  The  value  of  the  signal  am- 
plitude variance  had  significant  effect  upon  the  individual  station  capability 
but  little  effect  upon  the  total  network  detection  capability.  Finally,  an  eval- 
uation of  the  individual  stations  comprising  the  baseline  network  used  in  the 
study  yielded  a ranking  of  their  detection  capability. 

For  the  adaptive  beamforming  (ABF)  detector  study  using 
KSRS  short-period  data,  estimation  of  the  detection  probability  using  the 
beamsteer  and  ABF  detectors  was  accomplished  for  38  Eurasian  events  re- 
corded in  November  1974.  Preliminary  results  suggest  that  the  ABF  yields 
about  a 2 dB  improvement  over  beamsteering.  Currently,  a simulation  of 
burying  a scaled  signal  in  a given  noise  sample  is  being  undertaken.  This  is 
part  of  an  effort  to  improve  the  ABF  capabilities  for  weak  signals. 

A program  has  been  written  to  design  and  apply  a Wiener  filter 
to  long -period  signals  buried  in  noise  as  part  of  the  cascading  study.  Statis- 
tically reliable  noise  and  signal  power  spectra  have  been  found  in  the  course 
of  this  study,  and  the  signal-to-noise  ratio  gain  achieved  exceeds  that  of  a 
simple  bandpass  filter  by  as  much  as  3 dB  for  a realistic  case.  A more 
rapidly  adapting  version  of  the  three  component  adaptive  filter  has  also  been 
developed.  It  is  now  being  tested  on  several  events. 

During  the  last  quarter,  a programming  error  was  found  in  the 
upgraded  version  of  the  three  component  adaptive  processor  (TCA).  All  re- 
levant data  were  re-processed  to  correct  the  evaluation  of  the  TCA.  No  sig- 
nificant changes,  though,  were  found  in  the  evaluation  results. 
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The  dispersion  relation  filter  is  nearing  completion.  A time 
variant  Wiener  filter  algorithm  and  utilization  of  a variable  bandwidth  along 
the  dispersion  curve  have  been  incorporated  into  the  program.  This  Wiener 
filter  is  being  tested  on  synthetic  chirp  waveforms  and  actual  signals  buried 
in  seismic  noise.  The  filter  appears  to  work  reasonably  well  on  chirp  wave- 
forms. Also,  the  choice  of  the  expected  signal-to-noise  ratio  as  a Wiener 
filter  design  parameter  has  been  found  to  be  rather  critical. 

A Short-Period  Earthquake/ Explosion  Discriminator  (SPEED) 
processing  module  has  been  integrated  into  the  Interactive  Seismic  Process- 
ing System  (ISPS)  graphics  capability  on  the  PDP-15  computer.  A demon- 
stration of  SPEED  was  completed,  which  showed  capability  for  on-line  cep- 
strum  analysis  and  variable  frequency  magnitude  (VFM)  measurement.  These 
capabilities  were  demonstrated  on  a presumed  explosion  from  Kazakh  and  an 
earthquake  from  Iran  using  NORSAR  seismic  data.  The  results  strongly  sug- 
gest that  is  is  possible  to  carry  out  the  discrimination  processing  with  a sin- 
gle analyst,  keeping  up  with  the  anticipated  data  flow  of  25  to  50  events  per 
day  of  an  operating  surveillance  system.  The  SPEED  demonstration  showed 
that  cepstrum  analysis  could  be  done  to  reduce  false  identification  of  earth- 
quakes (as  explosions)  to  an  acceptably  low  rate  and  that  this  could  be  follow- 
ed by  VFM  processing  which  has  been  shown  to  have  about  as  much  discrimi- 
nation capability  as  Mg  versus  m^. 

Higher  mode  Love  waves  with  periods  between  5 and  10  seconds 
have  been  examined  during  the  last  quarter.  Studies  of  these  waves  generated 
from  two  Nevada  explosions  have  confirmed  that  the  energy  attenuation  coef- 
ficients for  the  source-station  (SDCS)  paths  are  relatively  high.  Because  of 
this  'anomalous'  attenuation,  accurate  measurements  must  be  made  in  order 
to  have  the  proper  correction.  Also,  the  three  component  adaptive  filter  was 
used  to  separate  higher  mode  arrivals  from  two  presumed  explosions  and  an 
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earthquake.  No  depth  discrimination  capability  was  found  by  fitting  these 
higher  mode  spectra  with  a theoretical  model. 
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SECTION  II 

ARRAY  AND  NETWORK  EVALUATION 

A.  ILPA  AND  SRO  EVALUATION 

1,  Current  Status 

During  the  past  quarter,  data  began  arriving  from  three  of  the 
Seismic  Research  Observatories  (SRO's).  These  reporting  SRO's  are  Albu- 
querque, New  Mexico,  Guam,  and  Mashhad,  Iran,  Processing  and  analysis 
of  long-period  and  short-period  data  from  these  SRO's  has  been  initiated, 
with  the  goal  of  establishing  the  noise  characteristics  and  signal  detection 
and  discrimination  capability  of  each  SRO.  From  the  limited  data  base  pro- 
cessed to  date,  several  preliminary  results  have  been  found. 

• Short-period  RMS  noise  as  recorded  at  Albuquerque  averages 
about  0.4  mfjL  (12  daily  samples). 

• Short-period  RMS  noise  as  recorded  at  Guam  averages  about 
30  m/x  , presumably  due  to  wave  activity  (12  daily  samples). 

• Long-period  RMS  noise  as  recorded  at  Albuquerque  averages 
about  9 m/u.  (14  daily  samples). 

The  following  preliminary  detection  capability  results  are 
based  on  a suite  of  events  with  epicentral  distances  predominantly  between 
60°  and  103°  for  short-period  data  and  between  60°  and  120  for  long-period 
data. 

• Considering  only  those  events  for  which  the  short-period  de- 
tector permitted  data  to  be  recorded,  the  short-period  50  per- 
cent detection  threshold  at  Albuquerque  is  at  approximately 

m = 4.  4. 
b 
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• Considering  all  events  in  the  data  base,  where  those  for  which 
data  was  not  recorded  are  counted  as  non-detections,  the  short- 
period  50  percent  detection  threshold  at  Albuquerque  is  at  ap- 
proximately m = 4.8. 

• Considering  all  events  in  the  data  base,  the  long-period  50 

percent  detection  threshold  at  Albuquerque  is  at  approximately 

m.  = 4.  6. 
b 

The  quality  of  the  SRO  data  sampled  has  been  in  general  excel- 
lent. On  some  samples  of  Albuquerque  data,  long-period  transients  (glitches) 
have  been  observed  on  the  north-south  component.  These  have  been  tenta- 
tively attributed  to  thermal  effects  by  personnel  at  the  Albuquerque  Seis- 
mological  Center.  No  spikes  have  been  observed  on  data  from  any  of  the 
three  SRO's  sampled  to  date. 

For  the  ILPA  evaluation,  the  Alaskan  Long-Period  Array  edit 
program  was  modified  to  adapt  it  to  the  Iranian  Long-Period  Array  satellite 
format.  The  output  of  this  modified  program  will  be  compatible  with  the  rest 
of  the  long-period  array  processing  package,  making  possible  a complete 
(noise  and  signal)  evaluation  of  ILPA. 

Testing  and  debugging  of  the  ILPA  edit  program  have  been  de- 
layed by  lack  of  test  data. 

Finally,  work  has  been  started  on  the  program  to  transfer 
ILPA  data  to  the  mass  storage  center, 

2.  Future  Plans 

During  the  next  quarter,  processing  and  analysis  of  both  long- 
period  and  short-period  data  recorded  at  the  previously  mentioned  SRO's 
will  continue,  building  up  a large  data  base  of  signal  and  noise  samples.  As 
other  SRO's  begin  to  report  data,  they  will  be  included  in  the  processing  and 

analysis  scheme. 


When  data  becomes  available,  the  ILPA  edit  program  will  be 
tested  and  debugged.  When  this  is  completed,  evaluation  of  the  noise  char- 
acteristics and  signal  d<  ection  capability  will  be  started. 

Also,  work  will  continue  on  the  ILPA  data  transfer  program. 
Much  of  the  time  will  be  spent  learning  the  data  language  necessary  for  the 
program. 

B.  AUTOMATIC  SIGNAL  DETECTOR  EVALUATION 

1.  Current  Status 

Major  accomplishments  for  work  on  the  optimum  automatic 
signal  detector  evaluation  include  the  following: 

• Detection  probabilities  using  37  Eurasian  events  recorded  in 
November  1974  were  obtained  for  several  prefilters  with  0.8 
and  1.  6 second  integration  gates.  Only  27  0°  and  300  azi- 
muths with  the  same  15  km/ sec  velocity  were  used  for  the  de- 
tector beams.  Figure  II- 1 shows  the  m^  distribution  (the  up- 
per part)  for  the  events,  and  the  maximum  likelihood  esti- 
mated detection  probability  curve  (the  lower  part)  for  the  0.  8 
second  integration  gate  Fisher  detector  result  (note  - OPTSS 
stands  for  optimum  detection  prefilter  designed  from  single- 
sensor data).  Complete  results  of  the  detection  probability 
from  the  event  ensemble  are  given  in  Table  II—  1 for  the  var- 
ious prefilters.  For  optimum  detection  with  the  Fisher  and 
the  conventional  power  detectors,  the  OPTSS  prefilter  has 
been  selected  for  further  study. 

• The  detector  response  pattern  was  studied  using  an  event  from 
Greece.  The  Fisher  and  conventional  power  detector  outputs 


DETECTION  PROBABILITY 


MB  DISTRIBUTION  OF  PROCESSED  EVENTS 
FILTER*  OPTSS  DETECTOR*  FISHER 


INTEGRATION  GATE* 


MAXIMUM  LIKELIHOOD  DETECTABILITY  CURVE 

FILTER*  0PTS5  DETECTOR*  FISHER  INTEGRATION  GATE' 


FIGURE  II-l 


THE  MAXIMUM  LIKELIHOOD  ESTIMATED  DETECTION 
PROBABILITY  CURVE  FOR  37  EURASIAN  EVENTS 
USING  THE  FISHER  DETECTOR 
II-4 


j] 

0 

0 


i 

I 

I 

i 


mm 

ektt 


i 

1 

I 

1 


os 

o 

Cm 

co 

w 

H 

►H 

►H 

CQ 

O 

OS 

CL 

z 

o 

HH 

H 

U 

W 

H 

_ W 
7 Q 
»— * 

~ Cm 

w o 
3 w 

< H 
H < 

S 

►H 

H 

co 

W 

Q 

O 

O 

S3 

►H 

X 

W 

K 

HH 

-I 

D 

2 


co 

3* 

<3 

OS 

3 

Q 
O 

•— i 

OS 

w 

CL  . . 

H OS 
OS  W 


r~ 


O 
S3 
co 

co 
OS 
co 
K 

H 
<3 

co 
OS 
O 
H 
U W 


PQ 

s 

w 

> 

o 

z 

z 

H-l 

co 

H 

Z 

W 

> 


£ 

<3 


3 

OS  o' 

WD 

> w 

O r~ 

CL  Pi 

SI  Cm 

2° 
£ co 

O ^ 

h3 
£ « 
W 


> 

£ 

o 

o 

Q 

£ 

< 

OS 

w 

S3 

co 

t-t 

Cm 


W 

S3 

H 

X 

o 


Ml 

o 

4-> 

') 

Mi 

4-> 

<L> 

Q 

Ml 

d 

* 

0 

Oh 

i-H 

d 

O 

•H 

4-J 

c 

d 

> 

d 

o 

U 


Ml 

o 

-Mi 

u 

a; 

-Mi 

d 

Q 

Ml 

d 

Si 

id 


i-H  i-H 

r-  tji 

vD  O 

O'  -H 

co  co 

ro  m 

i-H  i-H 

-H  00 

• • 

• • 

• • 

• • 

o o 

o o 

o o 

o o 

b 

+1  -H 

+!  -H 

+1  +1 

-H  +1 

o o 

o o 

00  ^ 

CO  Is- 

00  00 

00  00 

iji  m 

LD  LD 

o o 

0. 

0. 

o o 

o o 

00  00 

o m 

O'  v£> 

•-I  00 

ro  ^ 

CD  CD 

i-H  i-H 

oo  oo 

• • 

• • 

• • 

• • 

o 

o o 

o o 

o o 

o o 

O' 

JD 

-H  -H 

+1  +1 

-H  -H 

+1  +1 

c 

00  00 

00  o 

O O' 

O Tf 

G 

00  Tjl 

0-  o 

0-  \D 

o-  r- 

LD  LD 

4. 

5. 

TO 

TO  to’ 

O'  00 

00 

OO  o 

vO  t- 

H 

tn  oo 

i-H  i-H 

i-H  rt 

• • 

• • 

• • 

• • 

O 

o o 

o o 

o o 

o o 

LD 

_Q 

+1  +1 

-H  +1 

-H  -H 

-H  +1 

c 

vd  in 

m oo 

0s  in 

N ^ 

M 

00  ^ 

0-  0s 

o oo 

o o 

Tf  TO 

CO  CD 

Tj3  to 

to’  tO* 

i-H  sO 

O 0- 

00  o 

OO  O' 

ro  00 

OO  —I 

oo  00 

rH  00 

• « 

• • 

• • 

• • 

o o 

o o 

o o 

o o 

b 

-H  +1 

+1  M 

M +1 

+1  -H 

o m 

oo  in 

o o 

I-H  LD 

oo  r- 

so  in 

00  s> 

r- 

• • 

• • 

• • 

• • 

o o 

o o 

o o 

o o 

00  00 

00  00 

sO  m 

00  —1 

in  m 

oo  oo 

0O  00 

-H  CO 

• • 

• • 

• • 

• • 

o 

o o 

o o 

o o 

o o 

O' 

J3 

+1  M 

■H  +) 

+1  -H 

+1  +1 

g 

-h  in 

O' 

LD  LD 

00  00) 

00 

^ O' 

-H  O' 

sO  o 

• • 

• • 

• • 

• • 

LD  LD 

in  ^ 

in  to 

Tf  in 

O vO 

>-H  LD 

i-H  i-H 

DJ  ^ 

i-H  i-H 

00  -H 

-H  00 

• » 

• • 

• • 

* * 1 

o 

o o 

o o 

o o 

o o 

ID 

-O 

+1  +1 

+1  +| 

-H  -H 

+1+1 

6 

00  -H 

O Tf 

oo  r~ 

O 0- 

r-  LD 

00 

i-H  i-H 

— i o 

• • 

• • 

• • 

• • 

4 

4 

Tf  TO 

d 


0 

o 

•rH 

d 

nJ 

t/3 

u 

WD 

d 

<D 

HJ 

-Mi 

d 

d 

HH 

0 

Mi 

d 


d 

Mi 

X 


00  NO 

• • 

O — I 


00  sO 

• • 

O — I 


00  vO 

• • 

O -H 


N 

X 


3 


I 

ID 


X 

d 

d 

s> 

ID 

ID 

rt 

X 


II- 5 


d 

XI 


N 

X 

'fl  LO 
li  • 

■a  "7 
w ‘ 

tfi  CO 
n)  • 

X S 


6 


Ml  N 

d X 

d 

■ CQ  00 
d • 
^ •>,  00 

rt  ' 

M Mi  ® 

CM  Ml  • 

n < rL 


• • 

6 

d 

6 


00  sO 

• • 

O -H 


Mi  

0 TT 

03  S3 

..  d ® 

|»m. 

E J8  ^ 

•rt  OC  in 

a .5  j 

O co  ^ 


were  obtained  for  20  km/sec,  15  km/ sec,  and  10  km/ sec 
velocities  of  propagation  and  several  azimuths. 

The  response  pattern  was  examined  for  detector  output  loss  if 
a limited  number  of  beams  were  used  to  cover  a large  area  of  interest.  Pre- 
liminary results  suggest  that  the  Fisher  detector  yields  a larger  loss  than 
the  conventional  power  detector. 

2.  Future  Plans 

During  the  next  quarter,  detection  probabilities  for  three  re- 
gions will  be  examined  using  February  1976  (winter  season)  data.  Detection 
performance  using  single  sensor  data  will  also  be  studied. 


SECTION  III 

SIGNAL  DETECTION  METHODS 

A.  NEx  WORK  CAPABILITY 

1.  Current  Status 

A modified  version  of  che  program  NETWORTH  has  been  used 
to  evaluate  the  effect  of  certain  physical  parameters  on  network  detection 
capability.  This  study,  which  was  completed  during  the  last  quarter,  indi- 
cates that  the  station  detection  threshold  and  the  station  reliability  are  the 
primary  determinants  to  network  performance.  The  value  of  the  signal  am- 
plitude variance  had  significant  effect  ipon  the  individual  station  capability 
but  little  effect  upon  the  total  network  detection  capability;  a decrease  in  de- 
tection capability  at  one  station  was  generally  compensated  by  an  increase  at 
another  station.  Finally,  an  evaluation  of  the  individual  stations  comprising 
the  baseline  network  used  in  the  study  yielded  a ranking  of  their  detection 
capability. 

2.  Future  Plans 

The  report  on  this  evaluation  of  the  program  NETWORTH  has 
been  written  and  is  currently  being  reviewed  ,3rior  to  typing. 

B.  ADAPTIVE  BEAMFORMING  DETECTOR 

1„  Current  Status 

Estimation  of  the  detection  probability  for  38  Eurasian  events 
recorded  in  November  1974  using  the  beamsteer  and  ABF  detectors  were 
computed  with  the  maximum  likelihood  method.  Table  HI-1  shows  the  Gaussian 
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TABLE  III- 


m-2 


parameters  of  the  detection  probabilities  from  the  beamforming  and  ABF  re- 
sults. The  criteria  for  detection  was  that  the  signal  amplitude  was  at  least 
1.5  times  higher  than  the  noise.  These  preliminary  results  suggest  that  the 
ABF  yields  about  a 2 dB  improvement  ever  beamsteering. 

Also,  during  the  last  quarter,  the  array  beamforming  response 
pattern  was  examined.  Signal-to-noise  ratios  were  computed  by  beamforming 
an  event  from  Greece  with  20  km/ sec,  15  km/sec,  and  10  km/ sec  velocities 
for  various  azimuths.  The  results  of  this  analysis  will  be  discussed  in  the 
next  monthly  report. 

2.  Future  Plans 

Currently,  a simulation  of  burying  a scaled  signal  in  a given 
noise  sample  is  being  undertaken.  In  the  future,  detection  of  weak  signals 
will  be  examined  by  (1)  stopping  ABF  updating  upon  detection  by  the  Fisher 
detector,  and  by  (2)  using  low  ABF  convergence  rates  (~  0.  05). 
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SECTION  IV 

SIGNAL  ESTIMATION  TECHNIQUES 


A.  CASCADING  STUDIES 
1.  Current  Status 

A program  has  been  written  to  design  and  apply  a Wiener  fil- 
ter to  long-period  signals  buried  in  noise.  Statistically  reliable  noise  and 
signal  power  spectra  have  been  found  in  the  course  of  this  study,  and  the 
S/N  ratio  gain  achieved  exceeds  that  of  a simple  bandpass  filter  by  as  much 

as  3 dB  for  a realistic  case. 

A more  rapidly  adapting  version  of  the  three  component  adap- 
tive filter  has  been  developed.  It  is  now  being  tested  on  several  events. 


2.  Future  Plans 

The  effect  of  signal-noise  correlation  on  the  Wiener  filter  is 
going  to  be  investigated;  the  problem  is  modeling  this  correlation.  Also,  an 
averaged  matched  filter  waveform  will  be  found,  and  its  performance  and 
that  of  the  three  component  adaptive  filter  will  be  tested  alone  and  in  various 

combinations. 

B.  THE  THREE  COMPONENT  ADAPTIVE  PROCESSOR 
1„  Current  Status 

During  the  last  quarter,  a programming  error  was  found  in 
the  upgraded  version  of  the  three  component  adaptive  processor  (TCA).  All 
relevant  data  were  re-processed  to  correct  the  evaluation  of  the  TCA. 


IV- 1 


However,  no  significant  changes  in  the  evaluation  results  were  found.  The 
final  conclusions  of  the  evaluation  are: 

• The  most  effective  form  of  the  upgraded  Love  wave  TCA  pro- 
cessor contains  rotation  of  the  Love  wave  frequency  compo- 
nents about  the  vertical  axis  to  the  transverse  component  with 
an  accept-reject  limit  placed  on  the  Love  wave  arrival  azi- 
muth. 

• At  12  dB  true  signal-to-noise  ratio,  we  can  expect  7-9  dB  gain 
for  Rayleigh  waves  and  6-7  dB  gain  for  Love  waves  from  the 
original  TCA  processor  when  applied  to  single-site  data. 

• At  12  dB  true  signal-to-noise  ratio,  we  can  expect  about  9 dB 

* 

gain  for  Love  waves  from  the  upgraded  TCA  processor  when 
applied  to  single- site  data. 

• Only  low  or  negative  gains  can  be  expected  from  either  version 
of  the  TCA  processor  when  applied  to  beam  data. 

• Application  of  the  original  TCA  processor  to  single-site  data 

lowered  the  50  percent  Rayleigh  wave  detection  threshold  by 

about  0.25-0.40  m,  units  and  the  50  percent  Love  wave  detec- 

b 

tion  threshold  by  about  0.  15  m^  units. 

• Application  of  the  upgraded  TCA  processor  to  single-site  data 

lowered  the  50  percent  Love  wave  detection  threshold  by  about 

0.  35  m,  units, 
b 

• Application  of  the  original  TCA  processor  to  beam  data  lower- 
ed the  50  percent  Rayleigh  wave  detection  threshold  by  about 

0.  15  m,  units.  Essentially  no  change  occurred  in  the  50  per- 
b 

cent  Love  wave  detection  threshold. 
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• Application  of  the  upgraded  TCA  processor  to  beam  data 

yielded  essentially  no  change  in  the  50  percent  Love  wav^ 
detection  threshold. 

2.  Future  Plans 

The  report  on  the  evaluation  of  the  upgraded  TCA  processor 
has  been  written  and  is  currently  being  reviewed  prior  to  typing. 

C.  DISPERSION  RELATION  FILTER 

1.  Current  Status 

The  dispersion  relation  filter  is  nearing  completion.  A time 
variant  Wiener  filter  algorithm  and  utilization  of  a variable  bandwidth  along 
the  dispersion  curve  have  been  incorporated  into  the  program.  This  Wiener 
filter  is  being  tested  on  synthetic  chirp  waveforms  and  actual  signals  buried 
in  seismic  noise.  The  filter  appears  to  work  reasonably  well  on  chirp  wave- 
forms. The  choice  of  the  expected  signal-to-noise  ratio  as  a Wiener  filter 
design  parameter  is  rather  critical. 

The  filter,  to  some  extent,  is  capable  of  resolving  signal  esti- 
mates down  to  a -6  dB  signal-to-noise  ratio.  However,  it  can  also  generate 
false  signal  estimates  in  the  form  of  partial  chirp  waveforms  from  pure 
noise  which  are  difficult  to  discern  from  resolved  signals  up  to  +6  dB  signal- 
to-noise  ratio. 

2.  Future  Plans 

Testing  of  the  Wiener  filter  will  be  continued,  and  an  attempt 
will  be  made  to  solve  the  false  estimate  problem.  It  is  expected  that  the 
major  portion  of  the  technical  report  for  this  task  will  be  written  during  the 
next  quarter. 
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SECTION  V 
DISCRIMINATION 

A.  PDP- 15  DISCRIMINATION  PACKAGE 
1.  Current  Status 

A Short-Period  Earthquake/ Explosion  Discriminator  (SPEED) 
processing  module  has  been  integrated  into  the  Interactive  Seismic  Process- 
ing System  (ISPS)  graphics  capability  on  the  PDP- 15  computer.  A demon- 
. stration  of  SPEED  was  completed,  which  showed  capability  for  on-line  cep- 
strum  analysis  and  variable  frequency  magnitude  (VFM)  measurement.  These 
capabilities  were  demonstrated  on  a presumed  explosion  from  Kazakh  and  an 
earthquake  from  Iran  using  NORSAR  seismic  data.  The  results  strongly  sug- 
gest that  it  is  possible  to  c:.rry  out  the  discrimination  processing  with  a sin- 
gle analyst,  keeping  up  with  the  anticipated  data  flow  of  25  to  50  events  per 
day  of  an  operating  surveillance  system.  Allowing  a utilization  of  50%  of 
computer  capacity,  the  processing  time  per  event  for  all  routine  event  class- 
ification processing  should  be  on  the  order  of  15  minutes.  Specific  algorithms 
should  therefore  consume  no  more  than  several  minutes  per  event.  The 
SPEED  demonstration  showed  that  cepstrum  analysis  could  be  done  to  reduce 
false  identification  of  earthquakes  (as  explosions)  to  an  acceptably  low  rate 
and  that  this  could  be  followed  by  VFM  processing  which  has  been  shown  to 
have  about  as  much  discrimination  capability  as  Mg  versus  m^. 

One  of  the  major  problems  in  designing  a practical  interactive 
discrimination  processor  is  the  reduction  of  data  channels  to  one  or  several 
channels  per  event  (post  DAP  processing).  The  approach  followed  and  demon- 
strated was  to  apply  complex  cepstrum  beamforming  to  all  of  the  NORSAR 
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channels  to  produce  a single  event  channel  which  was  corrected  for  system 
response.  This  was  done  in  such  a way  that  the  analyst  could  choose  an 
event  channel  corrected  to  ground  displacement,  ground  acceleration,  or  the 
uncorrected  seismic  response.  Since  the  low  frequency  response  of  the  dis- 
placement wavelet  is  greatly  different  from  that  of  the  pP  operator,  the  best 
results  appear  to  be  obtained  by  choosing  the  channel  corrected  to  ground 
displacement. 

2.  Future  Plans 

SPEED  is  being  modified  to  include  an  absorption  correction  to 
the  spectrum  used  for  VFM,  Also,  the  deconvolution  program  which  separ- 
ates the  seismogram  into  direct  P and  echo  pP  wavelets  is  being  modified  for 
the  case  of  a filtered  echo.  Physically,  this  would  represent  a scattered  or 
multiple  echo.  Work  is  proceeding  to  make  the  post  DAP  data  reduction  an 
efficient  batch  process  suitable  for  production  usage.  A small  data  base  will 
be  processed  with  SPEED  to  evaluate  the  effectiveness  for  event  discrimina- 
tion. 


B.  HIGHER  MODE  STUDIES 
1.  Current  Status 

Higher  mode  Love  waves  with  periods  between  5 and  10  seconds 
have  been  examined  during  the  last  quarter.  Studies  of  these  waves  generat- 
ed from  two  Nevada  explosions  have  confirmed  that  the  energy  attenuation  co- 
efficients for  the  source -station  (SDCS)  paths  are  relatively  high.  Because  of 
this  'anomalous'  attenuation,  accurate  measurements  must  be  made  in  order 
to  have  the  proper  correction.  This  may  not  be  possible,  though,  with  the 
limited  number  of  stations  available. 

The  three  component  adaptive  filter  was  used  to  separate  high- 
er mode  arrivals  from  two  presumed  explosions  and  an  earthquake.  No  depth 
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discrimination  capability  was  found  by  fitting  these  higher  mode  spectra  with 
a theoretical  models  This  fitting  process  was  hampered  by  the  narrow  azi- 
muthal distribution  oi  the  SDCS  stations  in  relation  to  the  events  available 
for  analysis, 

2.  Future  Plans 

It  is  recommended  that  future  studies  for  this  task  involve  the 
possibility  of  studying  Eurasian  events  recorded  at  Eurasian  stations.  This 
should  reduce  the  problems  associated  with  attenuation  determination  and 
azimuthal  distribution  of  stations. 
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